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O I Abstract 
. 

Q^. We show that the second knee in the cosmic ray spectrum (i.e. the steep- 

^ \ ening occurring at -E ~ 4 x 10^^ eV) could be related to drift effects affecting 

\-{ ' the heaviest nuclear component, the iron group nuclei, in a scenario in which 

c/3 . the knee at 3 x 10^^ eV indicates the onset of drift effects in the lighter proton 

^ ' component. We also study the anisotropies resulting from diffusion and drift 

currents in the Galaxy, showing their potential relevance to account for the 
^ ' AGASA observations at E' ~ 10^^ eV, before the extragalactic component 

. becomes dominant. 

There are a couple of observations regarding cosmic rays (CRs) at energies ~ 
10^^ eV, i.e. just below the ankle, which still remain puzzling (see e.g. |l|, Q). One is 
the steepening in the spectrum taking place at an energy E^k — 4 x 10^'' eV, which 
is most apparent in the Yakutsk and Fly's Eye stereo 0, |^ data, and is usually 
referred to as the 'second knee', at which the spectral index a (such that dN/dE oc 
changes from a ~ 3 below Est to a ~ 3.3 above Est, although probably in a 
smooth way. The other is the 4% anisotropy at £^ ~ 0.8 — 2 x 10^^ eV reported by the 
AGASA collaboration , with an excess observed from a direction near the galactic 
center and a (probably associated) deficit near the galactic anticenter direction. Also 
some excess of events in the direction of the Cygnus region (along the Orion spiral 
arm) has been observed. An excess from a direction near the galactic center was also 
found in the SUGAR data but it has a smaller spatial extension (consistent with 
a point-like source) and its direction is few degrees away from the AGASA reported 
excess. The fact that these anisotropies are associated to galactic features is quite 
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important, because it is believed that the galactic cosmic ray component should just 
be fading away at > lO^^eV, since the CR spectrum should become dominated by 
the extragalactic component above the ankle, which lies at energies ~ 5 x 10^'^eV. 
Hence, the energies at which these features take place are probably the highest ones 
at which the CR spectrum is still dominated by the galactic component, and the 
features observed may actually be related with the very same process responsible 
for the fading away of this component. 

In this work we want to show that the steepening of the spectrum at Egk can 
be directly related to the enhancement of the escape mechanism of the iron group 
nuclei due to efficient drift effects, and that the observed anisotropics can be linked 
to the overall CR motion arising from these drifts. It is interesting that these 
processes would then be directly related to the scenario in which the knee in the 
spectrum at E^nee — 3 x 10^^ eV results just from the transition from a regime in 
which the CR transport is dominated by transverse diffusion to one dominated by 
the Hall diffusion, i.e. by drift effects ^. This scenario is based on the fact 
that each CR component of charge Z starts to be affected by drifts at an energy 
E c:=L Z Eknee, and its spectrum progressively steepens, with the spectral index finally 
changing by Aa ~ 2/3 in a decade of energy. The envelope of the total spectrum 
obtained by adding together the different nuclear components nicely fits the change 
from a spectrum oc E~'^''^ below the knee, to one oc i?"^ above it. However, above 
10^^ eV, where all the lighter components are strongly suppressed, the dominant 
iron component will progressively steepen its spectrum until the overall spectrum 
becomes oc E^^-^ above a few xlO^^ eV, hence reproducing the behavior observed 
at the second knee. This scenario also naturally explains the transition towards a 
heavier composition above the knee reported by several experiments. 

The change in the spectrum from the diffusion to the drift dominated regimes 
can be simply understood from the steady state diffusion equation V ■ J = Q, where 
Q is the source and the CR current is related to the CR density N through 

J = -Ds_W±N - DiiV\\N + DAhxVN, (1) 

with b being the unit vector in the direction of the regular magnetic field B^eg, i-e. 
b = Breg/|Breg|, and V|| = b(b- V), while = V — Vy. The components of the dif- 
fusion tensor are D\\ (along the direction of B^eg), D± in the perpendicular direction, 
while Da is associated to the antisymmetric part and determines the drift effects. As- 
suming for simplicity that the regular magnetic field is directed in the azimuthal di- 
rection and that in a first approach the Galaxy can be considered to have cylindrical 
symmetry, one finds that D\\ plays no role in the diffusion equation. Using now the 
fact that (for a Kolmogorov spectrum of fluctuations in the turbulent magnetic fleld 
component) D±{x.) ~ D\{-k){E / E^f'^ , while L'a(x) ~ D\{x.)E/Eo [|, one flnds 
that at low energies {E < Z Eknee) the transverse diffusion is the dominant process 
affecting the transport of CRs and leading to dN/dE oc {D^j_/ D±)dQ/dE oc 
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Figure 1: CR spectra computed for different galactic magnetic field models and 
source locations, which correspond to the same source and model parameters as 
in Figures 3 and 4 (see more details in the text). The dotted straight line is the 
extragalactic flux from ref. [1]. Also shown are the relevant experimental data points. 

where (3 is the spectral index of the source, while at high energies {E ^ Z Eknee) 
one has instead dN/dE oc {D\/DA)dQ/dE oc E'^-^ 

In Figure 1 we display the CR spectrum, computed following ref. for the 
galactic component, and adding an extragalactic component oc E~'^''^^ (with nor- 
malization taken from [|ll) to fit the observations beyond the ankle. Different curves 
represent different assumptions for the source distribution and for the galactic mag- 
netic field model. The adopted galactic source spectra are featureless extrapolations 
of the spectra measured below the knee, and hence the resulting changes in slope are 
just due to the energy dependence of the diffusion process responsible for the escape 
of CRs from the Galaxy. The choice of the model parameters was just guided by the 
requirement of correctly reproducing the first knee, and they are all quite plausible. 
As is apparent from the figure, the fit to the data is in general remarkably good. 

Also notice that if the extragalactic flux can be considered to be isotropic, it will 
not be enhanced by the Galactic diffusion process if no reacceleration in the Galaxy 
occurs. This may be understood as being a consequence of the Liouville theorem 



T0|| , and can also be seen from the diffusion equations by noting that the solution 
in the absence of sources and with isotropic boundary conditions at the outskirts of 
the Galaxy is just to have a constant density, so that there is no way in which an 



3 




Figure 2: Fractional abundances of different CR components, assuming that the 
extragalactic flux consists mainly of protons. The elements are grouped into light 
(1 < Z < 5), medium {Q < Z < 19) and heavy (20 < Z < 26) components. The 
data correspond to the same source and magnetic field model as in Figure 3. 

overdensity in the extragalactic flux could arise from diffusion processes alone. 

The relative abundances of different elements vs. energy are plotted in Figure 2, 
assuming that the extragalactic flux consists mainly of protons. One obtains that in 
our model the CR composition at -E ~ 5 x 10^'' eV consists of about 50% of galactic 
iron group nuclei (mostly Fe and Mn), of a proportion (~ 20%) of intermediate 
mass galactic nuclei [Z = 6-19) while the rest consists essentially of galactic and 
extragalactic protons. These numbers are quite consistent with the Haverah Park 
determination of nuclear abundances using the lateral distribution of air showers 
TT[| , which suggest that in the range 2 x 10^ ''-10^^ eV only ~ 30% of the CRs are 
light, with the rest being mostly heavy nuclear species. Also, the composition change 
to a lighter mix starting around the energy of the second knee is explained in 



this scenario as iron group nuclei drift away from the galactic plane with increasing 
efficiency for increasing energies. 

Turning now to consider the anisotropics resulting in this scenario, these will be 
produced by the galactic cosmic ray component alone, since as just discussed, under 
the assumption that at 10^*^ eV the extragalactic component is originally isotropic 
it will remain so after the propagation effects in the Galaxy are taken into account. 
The anisotropy associated to the galactic component is given by [1^ 

where is the total (galactic plus extragalactic) CR density, and J is the CR current 
introduced in Eq. (|l|). 
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Under the assumption of cylindrical symmetry and that the regular magnetic 
field is in the azimuthal direction, one has that V±N = WN and that the CR 
current will be perpendicular to the regular magnetic field (i.e. lying on the r — z 
plane). The contribution to the CR current arising from the transverse diffusion 
will be in the direction of VA^ (and hence orthogonal to the isodensity contours), 
while the drift part will be orthogonal to VA^ (and hence parallel to the isodensity 
contours). Under these simplifying assumptions, the components of the anisotropy 
will then read 



In Figures 3 and 4 we display, for different galactic magnetic field and source 
models, the density contours resulting from the numerical integration of the diffusion 
equations (solid lines) as well as the direction of the anisotropy vectors and with 
different shadings the amplitudes of the anisotropics. The results are for given E/Z 
values, i.e. for different energies according to which nuclear species are considered. 
Our location corresponds to tq = 8.5 kpc and zq ~ 0. The relative importance of 
the diffusion and drift components in the different galactic locations can be directly 
inferred from the relative directions between the isodensity contours and the arrows. 
Figure 3. a corresponds to energies well below the first knee {E/Z = 10^^ eV), where 
the diffusion is completely dominated by D±. The source is assumed to be in the 
galactic plane and to have a constant strength for r < 3 kpc, while vanishing at 
larger radii. We see that the resulting density contours are symmetric with respect 
to the galactic plane, and the anisotropics are perpendicular to them and quite 
small (~ 10~^). Figure 3.b is for an energy corresponding to E/Z ~ Eknee, where 
the drift effects start to be non-negligible and the local anisotropics are ~ 10"'^, while 
figure 3.C corresponds to rigidities such that E/Z ~ 1.5 x 10^^ eV/26, displaying 
the behavior that iron nuclei will follow for energies near 1.5 x 10^^ eV. In this 
regime the densities are largely determined by the drifts, and they refiect the general 
asymmetries of the regular magnetic field at large scales. The local anisotropics are 
in this case of order 10~^. 

Figure 3 was done for a regular magnetic field consisting of a disk component 
with reversals between arms, a vertical scale height Zd = 0.5 kpc and a local strength 
j^disk _ _Q 75 ^G0. The halo regular component had no reversals, a scale height 
Zh = 5 kpc, was symmetric with respect to the galactic plane and had a local value of 
j^haio _ _g_75 'jj^g random component had a similar scale height, Zr = 5 kpc 
and a local rms value of i?^""'^ = 1.5 /iG (see ref. for details of the magnetic 
field profiles adopted). In figure 4 we display results equivalent to those of figure 3.c 
but for different galactic magnetic field parameters and geometries and for different 
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Figure 3: Isodcnsity contours for different CR energies, with a) corresponding to 
E/Z = cV, b) to E/Z = Ekne,, and c) to E/Z = 1.5 x 10^^ cV/26. Every two 
contours there is a change by an order of magnitude in the densities. The arrows 
represent the direction of the anisotropy vectors, while the shadings indicate the 
amplitude of the anisotropics. The model corresponds to a source in the galactic 
plane constant inside 3 kpc, and the magnetic field structure described in the text. 

source locations r^. For instance, NR corresponds to a halo without reversals, while 
R (bottom left panel) corresponds to a halo with reversals similar to those adopted 
for the disk model. The label S{A) stands for symmetric (antisymmetric) halo 
models, i.e. such that Sj'^'°(r, z) = S^"'°(r, -z) {Bl^''^"{r, z) = -Sj"'°(r, -z)), while 
the + sign in the A+ label (top right panel) is because sign(i?Q"'° ■ B^^^^) is positive 
at our galactocentric radius for 2; > 0. Clearly the CR densities, the drifts and the 
corresponding anisotropics vary a lot from model to model, even if all models provide 
acceptable fits to the local CR spectrum, and hence it is seen that the observed local 
anisotropics can provide a useful handle to establish the general properties of the 
galactic magnetic field. 

It is important to notice that both the galactic center excess and the galactic 
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Figure 4: CR densities and drifts for different magnetic field models and source 
locations. In parenthesis are indicated the values of (Bq'^'', Eq"'", 5™"'^, za, Zh, z^) in 
//G and kpc respectively. 
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anticenter deficit observed by AGASA may be reproduced by an overall diffusion 
flow of the CRs near our location in the direction of the galactic anticenter, with 
an amplitude leading to a 4% anisotropy Since one should have reached at the 



relevant energies the regime in which the drift becomes dominant with respect to 
the transverse diffusion, and using the fact that locally siga^Bf^^) = —1, one should 
have as an approximate estimation that 

cN dz ' cN dr- ^ ' 

The overall local asymmetry will then point approximately in the radial direction if 
locally \dN/dz\ S> \dN/dr\. Its amplitude can be estimated from the expression of 
the antisymmetric diffusion coefficient, which is 

3 l + (c.r^P' 

where is the CR Larmor radius (with lo cjri^ the corresponding gyrofrequency) 
and the relevant time scale of velocity decorrelations (see ref. and references 
therein). This leads to 

I ^ {utaY ^ E 26kpc /iG {uta? 

' " 1 + (curA)2 " IQiSeVZ h \Breg\l + [ujtaY' ^ ' 

where h = \dln N/dz\~^ is the local vertical scale height of variation of the CR 
density. Let us notice that due to the mix of different nuclear species present, one 
has that actually 6 = J^ifi^i, with fi being the fractional abundances of galactic 
species of charge Zi (plotted in figure 2), while Si their corresponding local anisotropy 
at the energy under consideration. Although the presence of an extragalactic proton 
component tends to reduce the fractions fi, the nuclei lighter than Fe tend to have 
larger values of 6i (for some fixed energy), having then the tendency of increasing 
the value of 6. As a result, the values displayed in figures 3.c and 4, which are the 
anisotropics of Fe nuclei for 1.5 x 10^^ eV, are indeed good estimators of the total 
anisotropy 6 at this energy. 

We see from Eq. (^ that the anisotropics steadily grow with energy, reaching 
a few % at lO^^eV if the typical scale height of the CR density variations are of 
order a kpc. Significant vertical gradients in the CR density on the galactic plane 
naturally result in symmetric magnetic field models, since the vertical drifts in those 
models are similarly oriented in both hemispheres. To have the asymmetry vector 
pointing away from the galactic center clearly requires the local value of dN/dz to 
be negative, something which generally results when CRs have the tendency to drift 
to negative values of z, and this happens more pronouncedly e.g. in models where 
a symmetric halo field is pointed in the —0 direction. Also, in order that locally 
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\dN/dr\ < \dN/dz\ it is convenient that the source be not too close to our location. 
The model NR — A+ (top right) also shows the interesting feature that drifts tend 
to converge towards the galactic plane (the opposite would occur in A— models), 
resulting in an outward drift for z ^ 0. On the contrary, a halo model following 
the reversals of the disk (or a disk model alone without a regular halo), has very 
pronounced vertical drifts (bottom left panel). The factor {uta)'^/{1 + (utaY) in 
Eqs. (P) and (|^) reflects the suppression in Da in the presence of strong turbulence. 
It approaches unity if B^"'"''^ < B''^^ but it can become small if the turbulence is 
large, and this would then suppress the anisotropies. 

Let us mention here that other attempts have been made to explain the galactic 
center excess as resulting from CRs arriving almost rectilinearly from a central source 
|T^, |TB|, |T7|. Even if the CR were protons, the deflection produced by the regular (and 
random) galactic magnetic fields will largely exceed the 10° displacement between 
the observed excess and the galactic center direction (to produce this small deflection 
the source should be closer than ~ 2 kpc [|1§]). Hence, the hypothesis that the excess 
is due to neutrons produced from accelerated protons or nuclei through photopion 
production or spallation processes around the sources in the center of the galaxy has 
been explored as an alternative explanation P, [16|. The interesting fact about this 



proposal is that the decay length of neutrons is just ~ 10 kpc (£'/10^* eV), being of 
the order of the galactocentric distance for the energies relevant for the anisotropies, 
and hence this would also explain why no significant excess is seen at lower energies. 
Anyhow, to have such a powerful neutron source is not particularly naturalQ and 
furthermore there is no way in which such scenarios can also account for the deficit 
observed in the anticenter direction. 

As regards the excess of events observed from the Cygnus direction, let us notice 
that a similar explanation to that discussed here can in principle be searched in 
terms of an overall CR diffusion along the spiral arm. This would clearly require 
a departure from the assumption of cylindrical symmetry and to take into account 
the parallel diffusion of CRs along the arms. A general prediction of this kind of 
scenario would be that a deficit associated to the opposite direction would result, 
something of potential interest for future observations in the southern hemisphere 
with the AUGER observatory. 
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^On the contrary, the maximum energy of galactic protons in our scenario needs not be larger 
that ~ 10"'^'^ eV, a value which is clearly much more plausible. 
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